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Введение. Рассматриваются ракета-носитель (РН) в полете и динамические составляющие нагрузок от 
воздействия порыва ветра трапецеидального профиля. Предлагается определить динамические составляющие 
силовых факторов с помощью аналитических решений для ускорений точек конструкции. Цель работы — 
создание методики выбора продолжительности нормативного порыва, при воздействии которого достигаются 
максимальные нагрузки в сечениях корпуса ракеты-носителя. 

Материалы и методы. Ракета-носитель представляется в виде неравномерной балки. Описание ее колебаний 
приводится к системе независимых обыкновенных дифференциальных уравнений, определяющих движение 
эквивалентной системы осцилляторов. Уравнение колебаний осцилляторов под действием трапецеидальной 
импульсной нагрузки решается методом наложений и сводится к вычислению интеграла Дюамеля. 
Предлагается получать параметры эквивалентной системы осцилляторов по результатам расчета динамических 
характеристик для конечноэлементной модели РН в программе Мазгап. 

Результаты исследования. Приведены аналитические соотношения для ускорений точек корпуса РН под 
воздействием порыва ветра трапецеидального профиля. Для балочной модели проведены тестовые расчеты 
ускорений по методике, предложенной в настоящей статье. Эти данные сравниваются с результатами 
конечноэлементного моделирования. С помощью аналитических решений построены зависимости, 
определяющие характер изменения величины изгибающего момента для различных сечений ракеты-носителя 
при варьировании продолжительности порыва ветра. 

Обсуждение и заключения. Благодаря представленной методике можно на примере ракеты-носителя строить 
эквивалентную динамическую модель систем с большим числом степеней свободы и получать аналитические 
решения для ускорений точек механической системы при трапецеидальном внешнем воздействии. Данные 
решения применимы для исследования динамических нагрузок. Итоги такого анализа позволяют выбрать 
продолжительность порыва ветра, при которой достигаются максимальные нагрузки в сечениях корпуса РН. 
Вычисления на базе аналитических решений весьма экономичны с точки зрения затраченного времени. Их 
можно использовать в проектных расчетах для предварительной оценки нагружения. 
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Введение. При эксплуатации ракеты-носителя (РН) в элементах ее конструкции возникают нагрузки. 
Речь идет о продольных и перерезывающих силах, изгибающих и крутящих моментах. Данные об указанных 
силовых факторах задействуют для прочностного анализа при проектировании новых изделий, 
экспериментальной отработке конструкции [1] и адаптации средств выведения под конкретный пуск [2]. 
Нагрузки подразделяются на квазистатические и динамические. Квазистатические, возникающие в полете, 
вычисляются из условия динамического равновесия РН как твердого тела с учетом допустимых параметров 
траектории выведения. 

Для расчета динамических нагрузок обычно применяются обобщенные балочные модели. Такие 
нагрузки определяются по результатам решения уравнения движения упругой РН, которое в общем случае 
представляет собой уравнение в частных производных. Методы, основанные на использовании разложений 
решений по тонам собственных колебаний конструкции, способны обеспечить высокую скорость вычисления в 
сочетании с достаточной точностью результатов [3]. Как показано в [4], с помощью метода разложения по 
формам собственных колебаний можно перейти к системе независимых обыкновенных дифференциальных 
уравнений. Они описывают: 

— движение центра масс РН, 

— поворот продольной оси РН относительно центра масс, 

— поперечные упругие колебания РН в полете. 

Переход к независимым уравнениям, описывающим упругие колебания РН, означает, что будут 
учитываться распределенные параметры осцилляторов, каждый из которых представляет собой систему с 
одной степенью свободы. Движение такого осциллятора можно рассматривать независимо от других и для 
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каждого из них можно получить решение с помощью известных методов теории колебаний. 

С целью расчета нагрузок в полете рассматривается не вся траектория выведения, а только некоторые 
ее точки, так называемые случаи нагружения, характеризующиеся экстремальным значением отдельных 
параметров, влияющих на нагружение, либо максимальным значением нагрузок на отдельные элементы 
конструкции. Один из наиболее важных случаев нагружения — полет РН в плотных слоях атмосферы [1,5]. 
Влияние неспокойной атмосферы на нагружение РН можно определить статистическими методами [6, 7] либо в 
рамках консервативного подхода, когда учитываются максимально возможные (с некоторым уровнем 
вероятности) характеристики ветра. В данной статье рассматривается второй подход. В качестве внешнего 
динамического воздействия принимается однократный нормативный порыв ветра. Профиль нормативного 
порыва, характеризующий изменение скорости ветра с течением времени, может быть задан в 
трапецеидальном [8], косинусоидальном или синусоидальном виде [9]. В данной статье рассмотрим движение 
ракеты-носителя под действием трапецеидального порыва. Продолжительность нормативного порыва обычно 
выбирается сопоставимой с периодом низшего поперечного тона колебаний РН'. При этом часто выдвигаются 
требования ее варьирования для достижения максимальных усилий в сечениях РН [8, 10]. Трудоемкость 
расчетов с использованием стандартных программ конечноэлементного (КЭ) анализа обусловлена: 

— необходимостью варьирования параметров внешнего воздействия, 

— большим числом расчетных случаев, 

— множеством вариантов исполнения и конфигураций конструкций на этапе проектных расчетов [9]. 

Цель настоящей статьи — разработка методики выбора продолжительности действия нормативного 
порыва с использованием аналитических решений, полученных для упрощенной динамической модели РН, 
представленной в виде эквивалентной системы осцилляторов. 

Полуаналитический подход с использованием интеграла Дюамеля был успешно применен в [11] для 
проведения гидроупругого анализа судов. В [12] и ряде других работ интеграл Дюамеля используется в рамках 
решения задачи нагружения мостов подвижными нагрузками. В настоящей статье интеграл Дюамеля 
задействовали для аналитических решений реакции ракеты-носителя на кратковременное воздействие порыва 
ветра в полете. 

Материалы и методы. На этапах эскизного проектирования целесообразно использовать плоские 
расчетные схемы для балочных моделей. При простоте и скорости решения они позволяют определить 
параметры движения и внутренние усилия (с приемлемой для этого этапа проектирования точностью) [13]. 
Представим РН в виде упругой балки с переменными по длине массой и жесткостью. Примем обычные для 
сопротивления материалов допущения, в том числе гипотезу о малости упругих деформаций. Для определения 
внутренних силовых факторов в сечении РН используем метод ускорений (перегрузок) [1, 4], который можно 
интерпретировать как метод сечений, адаптированный для динамического расчета. В этом случае внутренние 
усилия находят из условий статического равновесия мысленно отсеченных частей конструкции под действием 
внешних распределенных нагрузок, дополненных силами инерции Д’Аламбера, и искомых внутренних усилий. 
Квазистатические и динамические значения силовых факторов определяются отдельно на основе 
предварительно вычисленных ускорений, а затем суммируются [4]. 

В данной работе рассматривается вопрос определения динамического нагружения РН в поперечном 
направлении под действием порыва ветра, скорость которого направлена перпендикулярно продольной оси РН. 
Принимается, что нагружение в продольном направлении можно рассчитать независимо. В данной работе оно 
не рассматривается. 

С целью определения динамических ускорений РН представляется в виде свободной упругой балки. Ее 
движение изучается в окрестности момента времени, соответствующего рассматриваемому случаю нагружения, 
и описывается в отклонениях от состояния динамического равновесия, в котором находилась РН до ветрового 
порыва, двигаясь по номинальной (невозмущенной) траектории. При этом как бы «замораживаются», 
т.е. принимаются постоянными и равными характеристикам рассматриваемой точки номинальной траектории 
такие параметры, как масса и момент инерции РН, тяга двигателя, угол наклона траектории. Возмущенное 
движение упругой РН исследуется в неподвижной системе координат, связанной с положением, которое РН 
занимала на момент начала расчета. Возмущенное движение будет представлять собой совокупность 
плоскопараллельного движения РН как твердого тела в плоскости, в которой приложена динамическая 
нагрузка, и упругих движений корпуса. Реакцию системы управления не учитываем, т.е. для автомата 
стабилизации предполагаем большое время запаздывания по сравнению со временем приложения 
динамической нагрузки. Вообще, под воздействием порыва ветра вместе с упругими колебаниями корпус РН 
начинает движение как твердое тело. Поперечная составляющая аэродинамической силы, которая при малом 
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угле атаки считается ему пропорциональной, меняет значение за счет смещения корпуса в направлении 
действия порыва ветра и поворота относительно вектора набегающего потока. Проекция силы тяжести на 
поперечную ось РН также меняется. Для правильного учета указанных изменений уравнения движения РН 
необходимо интегрировать с уравнениями, описывающими логику работы автомата стабилизации, что 
невозможно на ранних этапах проектирования. Учитывая к тому же значительную массу и момент инерции РН, 
будем считать малыми: 

— угол поворота РН как твердого тела за время расчета, 

— скорость смещения центра масс РН в направлении действия порыва ветра. 

Это позволяет не учитывать влияние указанных выше изменений. Приращение угла атаки (а 
следовательно, и поперечной составляющей аэродинамической силы) считается зависящим только от величины 
скорости порыва ветра, заданной в виде функции времени. Таким образом, с учетом принятых допущений 
аэродинамическая нагрузка в поперечном направлении представляет собой распределенную по длине балки 
нагрузку с коэффициентом пропорциональности, зависящим от времени. Закон распределения 
аэродинамической нагрузки по длине РН определяется экспериментально и считается заранее известным. Закон 
изменения коэффициента пропорциональности (угла атаки) от времени определяется выбором профиля 
нормативного порыва ветра. 

Движение РН, смоделированной в виде упругой балки, можно описать, используя известное уравнение 
вынужденных поперечных колебаний балки, записанное с учетом гипотезы Фойгта: 


то В + (1412) В] = (к, 0), (1) 


где т(х) — погонная масса; В(х) —  изгибеля жесткость; 4(х,#Ё) — распределенная внешняя нагрузка; й — 
коэффициент трения. 

Данное уравнение должно быть дополнено граничным условием: внутренние усилия в начальном и 
конечном сечении равны нулю. Значит: 


> | В( ее = 0, Во = 0, прих = 0,х = [, (2) 


где [ — длина РН. 
В качестве внешней распределенной нагрузки в данной работе принимается аэродинамическая сила, 
которую можно представить в виде произведения функций: 

4(х, 2) = В (У. (х), (3) 
где В (Е) — функция, определяющая временную изменчивость аэродинамической силы и изменяющаяся по 
трапецеидальному закону в соответствии с принятой в настоящей работе моделью порыва ветра; У„(х) — 
функция распределения аэродинамической силы по длине РН. 

Рассмотрим свободные колебания РН без учета сил трения (при 4(х, Е) = 0, Й = 0). Подставим метод 
разделения переменных У(х, Е) = Г(х) :9(Е). В этом случае от уравнения (1) с граничными условиями (2) 
можно прийти к обыкновенному ыы уравнению 


[Во] — рат) = 0, (4) 





ах? 


с краевыми условиями: 
Ч [В Г] = 0, В(х 2-е? = 0, прих = 0,х = [. (5) 


Решение (4) с условиями (5) представляет собой классическую задачу Штурма — Лиувилля. Решая ее, 
можно найти набор собственных форм Д(х) и собственных частот р; рассматриваемой балки (} = 1,2,...). 
Известно”, что некоторые решения системы (4) соответствуют нулевым собственным частотам. Формы, 
соответствующие нулевым собственным частотам, определяют поступательное движение РН как твердого тела 
вместе с центром масс и вращение вокруг центра масс: [1 = 1, К =Х - хс, где хе — координата центра масс 
РЕВ 

Следует отметить, что массовые и жесткостные характеристики РН чаще всего имеют кусочно- 
постоянный характер распределения. В этом случае уравнения вида (1) и (4) должны быть записаны отдельно 
для каждого однородного балочного участка с граничными условиями на стыках участков, как при выводе 
соотношений метода начальных параметров [1, 14]. В данной работе указанная запись опущена, т. к. расчет 
динамических характеристик (модальный анализ) конструкции проводится численно, с использованием метода 
конечных элементов. 

Представим вынужденные колебания упругой балки, моделирующей корпус РН, в виде разложения по 
формам собственных колебаний. Примем, что ось жесткости балки проходит через ее центр масс. Для 


? Колесников К. С. Динамика ракет. М. : Машиностроение, 2003. 520 с. 
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перемещения точек оси РН запишем: 

у(х, В = ус(® + 9(0(х- хо +. 0 а; (0, (6) 
где ус(Е) — перемещение центра масс балки; 9(Е) — угол поворота оси недеформированной балки; /; — форма 
собственных колебаний балки, соответствующая тону с номером ]; 4;(®) — обобщенная координата, 


соответствующая тону с номером /; М — число учитываемых упругих тонов. 
После подстановки (6) в (1) и применения процедуры Бубнова — Галеркина можно прийти к системе 
обыкновенных дифференциальных уравнений с постоянными коэффициентами: 


туУс(® = О0,, 
19(0 = 0... (7) 
т;(4;(® + 2п,4;(®) + рта) О) 
Здесь т — масса РН; Г — момент инерции относительно оси, проходящей через центр масс РН 
перпендикулярно плоскости вращения; т; — приведенная (обобщенная) масса для -го тона колебаний и 


определяемая формулой: т; = | в т(х) {7 о) ах. 


Обобщенные силы в выражении (7) с учетом (3) определяются следующим образом: 


р 
Чу = В (Е) р о ах = ВАО; (8) 
р 
0 = В@) ]› У, (®) (х — хсдах = В(®Ма, (9) 
р 
0; = В(Е) И Уз (х) рах = В(®)0у,. (10) 
Здесь О. — максимальное значение главного вектора поперечной аэродинамической нагрузки; Ма — 


максимальное значение главного момента поперечной аэродинамической нагрузки, приведенного к центру масс 
РН; Оо; — максимальное значение обобщенной силы, соответствующей обобщенной координате 4. 

Первые два уравнения в (7) определяют закон изменения ускорений точек РН в процессе 
поступательного и вращательного движения РН как твердого тела. Последнее уравнение в (7) определяет закон 
движения эквивалентной системы осцилляторов. 

Рассмотрим движение одного осциллятора под действием трапецеидальной внешней нагрузки, которая: 

— возрастает от нуля до (0% за время 6, 

— сохраняет постоянное значение за время 6, 

— падает до нуля за время 6. 

Для удобства опустим индексы, характеризующие номер тона. Тогда дифференциальное уравнение 
движения осциллятора с учетом (10) перепишем в виде: 


4+ 2п4 + ра =3 = (В), (1 
т т 

где т — масса осциллятора; р — угловая частота собственных колебаний, выраженная в радианах в секунду; и 

— коэффициент затухания (определяет демпфирующие осциллятора); О(Е) = О,В(Е) — закон изменения 

внешней нагрузки. 


Представим функцию В (Е) как совокупность четырех линейных функций: 


В (Е) =6 В2 (©) =-6 "(Е - 6), 
Вз(Р) =-6 "(&- 6-0), В. (Е) =6 "(Е- 26-0). 
Соответственно, внешняя нагрузка О(Г) — это совокупность четырех линейных нагрузок О;(Е) = 


ОВ; (Е), (Е = 1,2,3,4). Нагрузка 01 (Е) прикладывается с момента & = 0; 05(Е) — с момента & = 6; О03(Е) — с 
момента & = 0+5; О0.,(Е) — с момента Е = 0 + 265. 

Для определения реакции системы на внешнее воздействие разобьем все время действия нагрузки на 
четыре интервала (рис. 1). 


9(е),Н 





‚1-й участок ! 2-й участок '3-й участок! 4-й участок 


Рис. 1. Схема внешней трапецеидальной нагрузки 


Малыхина О. И. Аналитическое решение приближенных уравнений движения ракеты-носителя под действием порыва ветра 





В соответствии с методом наложений? 


представим в виде интеграла Дюамеля реакцию 
рассматриваемой линейной системы на внешнее воздействие как сумму реакций на совокупность независимо 


приложенных элементарных импульсов: 
Ё 


а(® = | Р(Е — 9)У(9)а8, 


0 
здесь Р(Е — 9) — смещенный временной закон изменения внешнего воздействия, а У (9) характеризует реакцию 


системы на единичное импульсное воздействие. Реакцию механической системы 4а1(Р) на линейно 
возрастающую нагрузку Р(Е) = КЕ можно выразить через реакцию системы на единичный импульс У(9) и на 


внезапно приложенную единичную нагрузку У' (9): 
Е 


41(Е) = К | у (9)а9, У(®= | У(9)а9. 
0 


0 
Реакция на единичный импульс для механической системы с одной степенью свободы, имеющей 


коэффициент затухания п, будет иметь вид": 
- И 
У (9) = ———е а ( р? — п28. 
Введем обозначение для частоты затухающих колебаний р. = \/р? — п?. Вычислим реакцию системы 


на внезапно приложенную единичную нагрузку: 
Е 


у. (9) = | (3) аз = 


0 
Реакция на линейно возрастающую нагрузку: 


41 (Е) = К Л у, (9)а9 = — Е [-1 +е-"* со$(р.®)] + (= — 1) е-"* зт(р1 Е) + Е. (12) 


тр? р 
Примем обозначения: Д.= # — 6, Д.= Д. — 0, Д.=Д. — 6. 





—- Е —е т (созфь5 -- 5-5 (210) | 





Введем функцию, содержащую гармонические члены решения (12): 
п п 1 
4(9) = 2— (собр .9) + — зт(ф: . 5) е-"® — — зт(р, 9) е"%. 
р Р1 Р1 


Обозначим 4‹. перемещение под действием статически приложенной к системе силы Оу = а-тр? и 
учтем, что К = 05/6. Суммарное перемещение осциллятора под действием комбинации нагрузок О;(Е) на 
каждом временном участке будет суммой соответствующих решений (12). Реакция системы, описываемой 
уравнением (11), на внешнее воздействие трапецеидального профиля будет иметь вид: 


Чстб а(® — 21/° +1 при0 < ЕЁ < 6, 
и тб [А(®) - а(4,) +6] приб < < 0+6, 3 
«6 Ча(® - а(А.) — а(4,) + 2п/р? —-4.|] приб +6 <Е< 0-25, 
Чсб`"(4(®) — а(№,) — а(4») + а(Аз)) пра 
Закон изменения ускорений осциллятора может быть получен двойным дифференцированием по 
времени выражения (13). Введем функцию 9(9) = “ (а(9)). Дифференцируем и запишем результат с учетом 


номера тона колебаний: 


п, п: 1 
9;:(9) = © — р?;) (7 [совр .9) + — (ри ; :9)) — —$т(р! ; :)) — 
р; Р1 у Р1 у 


АЙ, (27а, зт(ра ; .9) + п; соз(р: .9)) == соз(р: .9))) Е". (14) 


Перепишем выражение для частоты затухающих собственных колебаний: 


ра; = № Я. (15) 


Для определения ускорений точек РН дважды продифференцируем (6): 
У(Ех) = Ус (® + 9(@)(х-хо +. (ао. (16) 
Выразим ускорения обобщенных координат из первых двух уравнений системы (7) с учетом (8) и (9). 
Представим ускорения точек оси РН, которые она приобретает, двигаясь как твердое тело: 


а(®) = ус (® + 9(0(х — хе) = В((Оа/т + Ма/П. (17) 


°Бидерман В. Л. Теория механических колебаний. М. : 0ВЗ$, 2017. 416 с. 
Яблонский А. А. Курс теоретической механики. М. : Интеграл-Пресс. 2007. 603 с. 
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Кроме того, учтем закон изменения функции В (Е) от времени: 


5-1 при 0 < Е < 6, 
в(® = 1 приб < < 90 +65, 
—6 1(- 26-0) при +6 << 0+ 25, 
0 при & > 09 + 25. 
Получим: 
9 1 (Оа/т + М./П при 0 < ЕЁ < 5, 
а(® = От приб < & < 9+6, (18) 
—5-1(&-— 26-0)(0./т-+ М./Г) при9+б<Е<О+ 256, 
0 при > 0 + 26. 


Для ускорений обобщенных координат, соответствующих упругим тонам колебаний, в результате 
двойного дифференцирования выражения (13) с учетом (14) получим: 
6 Ч ;9 (Е) при О=Е< 6, 


6 "Чет; (я® — 9.) приб < < 0+6, 
5-14ь; (9/0 — 9.) - 9, 45)) прив +8<Е<0+25, 


5-14; (9/(® — 9, (В) — 9,4») + 9,(43)} приё>0+25. 

С учетом (16) и (17) закон изменения ускорений точек упругой оси РН под действием порыва ветра 
трапецеидального профиля будет иметь вид: 

УЕ, х) = а(® +=, 40, (20) 
где а(Е) определяется выражением (18), 4;(Е) — выражением (19). 

При исследовании упругих колебаний РН демпфирование традиционно принимается на основании 
данных, полученных по результатам натурных динамических испытаний и представленных в виде значений 
логарифмических декрементов Д;. Тогда коэффициент, определяющий параметр демпфирования и входящий в 
выражения (14) и (15), можно вычислить по формуле: 

п; = р г (21) 

Входящее в формулу (19) изменение обобщенной координаты 4; под действием статически 


4; (2) = (19) 


приложенной обобщенной силы Оо; определяется выражением: 
1 901 _ а Ха г, ах 
р} т — ру ут?) ах. 
Принимая во внимание (14), (15), (18), (19), (21), (22), формула (20) — это аналитическое выражение, 
определяющее функции изменения ускорений по времени для точек оси РН под действием внешней 
поперечной аэродинамической силы, изменяющейся по трапецеидальному закону с учетом влияния 
диссипативных сил. 
Зная закон изменения ускорений, можно по известным методикам [4] определить динамические, а 
затем и суммарные корпусные нагрузки, действующие в сечениях РН. Вычислим изгибающий момент, 
обусловленный инерционными силами, возникающими за счет упругих колебаний корпуса РН: 


Чст  — 22} 


М(х, в =-»/-1М(х)4, (0. (23) 

Здесь М;(х) — функция распределения по длине РН единичного (при ускорении 4,;(#), равном единице) 
изгибающего момента для -го тона колебаний. Ее можно найти по формуле: 

М; (х) = Г. Г то родах Ах. (24) 


Как отмечалось выше, для проведения модального анализа в данной работе использовался метод КЭ. 
Данный подход обусловлен тем, что на практике у динамической модели РН достаточно сложная структура. 
Она включает подконструкции, и нельзя пренебрегать их собственной динамикой. Подконструкции могут 
крепиться к корпусу РН в одном сечении либо располагаться параллельно продольной оси РН и иметь 
несколько точек крепления. В этом случае расчет динамических характеристик в континуальной постановке 
является сложной математической задачей. Кроме того, динамические модели отдельных подконструкций 
представляются предприятиями-разработчиками в конденсированном (матричном) виде в формате Мазгап. По 
этой причине для расчета динамических характеристик конструкции будет оптимальным задействовать пакет 
программ инженерного анализа Мазгап. 

Однако применение стандартных программ КЭ-анализа для расчета динамического нагружения 
сопряжено с определенными трудностями. К ним относится необходимость предварительно строить 
эквивалентную модель внешней аэродинамической нагрузки, пригодную для использования в программе КЭ- 
анализа [15], что само по себе достаточно сложно. Кроме того, есть трудности, связанные с обработкой 
результатов расчета. Применение для анализа результатов функционала постпроцессора чрезвычайно 
трудоемко, требует большого числа ручных операций. Другой путь предполагает применение дополнительного 
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программного обеспечения для обработки массивов числовых данных большой размерности”. В данной работе 
предлагается подход, в рамках которого программа КЭ-анализа используется только для модального анализа. 
Динамическое нагружение в этом случае рассчитывается с помощью специально разработанного программного 
обеспечения, позволяющего варьировать внешние нагрузки и автоматически обрабатывать результаты расчета. 
Для сложной конструкции РН уравнения (7) сохраняют свой вид [3]. Стандартная информация вывода 
программы Мазгап может быть базой для получения параметров эквивалентной системы осцилляторов, массы 
и момента инерции РН, а также для вычисления обобщенных сил, входящих в третье уравнение системы (7). 
Для формирования левой части третьего уравнения в (7) необходимы значения собственных частот р; (КаЧ1апз) 


и обобщенных масс т; (Сепега|2е4 таз). Для определения обобщенных сил в правой части третьего 
уравнения в (7) необходимы функции собственных форм колебаний }.(х) (Евепуесюг). При вычислении 


динамических инерционных нагрузок по (23) вместо (24) удобнее воспользоваться единичными инерционными 
нагрузками (силами и моментами), которые выводятся программой Мазгап после стандартного запроса усилий 
при расчете собственных форм и частот. Единичные инерционные нагрузки выводятся отдельно для каждого 
тона колебаний, причем помноженные на квадрат собственной частоты, что следует учитывать для их 
корректного использования. 

Результаты исследования. Для проведения тестовых расчетов рассмотрена ракета-носитель среднего 
класса тандемной схемы. Вид функций погонной изгибной жесткости В(х) и массы т(х), а также 
распределения сосредоточенных масс Теоср = Тсоср,А(х —х.) по длине рассматриваемой РН представлен 


на рис. 2 (через Д обозначена дельта-функция Дирака). 


т(х), кг/м 


Тсоср (х), кг 


ВН | | х, м 





в) 


Рис. 2. Распределение массовых и жесткостных характеристик по длине РН: 
погонная масса (а), сосредоточенная масса (6), изгибная жесткость (6) 


° Малыхина О.И. Автоматизация обработки результатов конечноэлементного анализа нагружения конструкций ракетно-космической 
техники : сб. ст. УП науч.-тех. конф. молодых ученых и специалистов Центра управления полетами. Королев : ЦНИИМаш, 2017. С. 427-— 
434. 
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Для проверки полученных аналитических решений при воздействии трапецеидальной внешней 
нагрузки рассчитаны динамические ускорения РН по формуле (20) с учетом (14), (15), (18), (19), (21), (22). При 
этом учтены 5 упругих тонов собственных поперечных колебаний РН. Характеристики р,;, т; и }/(х) были 
получены по результатам расчета в программном комплексе МЗС Мазбап с использованием 
последовательности решения для модального анализа собственных колебаний (ЗОГ. 103). Динамическая 
конечно-элементная модель РН представлена в виде набора балочных элементов с различными инерционными 
и жесткостными характеристиками. К балочным упруго или жестко присоединены: 

— элементы, описывающие инерционные свойства приборов, агрегатов, частей конструкции блоков; 

— конденсированные модели отдельных блоков, представленных в цифровом матричном виде. 

Конечноэлементная модель РН включает около 1000 составляющих. Элементы, моделирующие 
закрепления, не использовались для сохранения способности движения РН как твердого тела. 

По той же конечноэлементной модели РН в программном комплексе МС Мазиап рассчитаны 
ускорения с использованием последовательности решения ЗОГ, 119, применяемой для модального анализа 
переходных процессов. При этом в модальном разложении были учтены все тона собственных колебаний в 
диапазоне до 100 Гц. Аэродинамическая нагрузка представлена поперечными погонными нагрузками, 
распределенными по всем балочным элементам, моделирующим корпус РН. Кроме того, задан 
трапецеидальный закон изменения аэродинамической нагрузки по времени. 

На рис. 3 представлены результаты сравнительного анализа полученных с помощью двух различных 
подходов ускорений ао(Е) = У(Е, ху) некоторой точки оси РН с координатой х = хо. Видно, что хорошо 
согласуются два решения: 

—щ полученное для ускорений по упрощенной модели РН на основе аналитических соотношений, 
приведенных в настоящей работе; 

— численное, полученное по полной конечноэлементной модели РН. 
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Рис. 3. Сравнение поперечных ускорений, полученных с помощью двух различных подходов: — конечноэлементное 


моделирование, —_ аналитическое решение для упрощенной модели 


Выбранное из временного процесса пиковое значение изгибающего момента определяет уровень 
эквивалентных усилий, принимаемых для проведения прочностного расчета, и выступает изменяющимся 
параметром при варьировании параметров внешнего воздействия в поперечном направлении [16]. 

На основе значений динамических ускорений получены зависимости (рис. 4). Они показывают, как 
изменение пиковых значений изгибающего момента М для различных сечений РН связано с параметром 0, 
характеризующим продолжительность действия порыва. На рис.4 значения изгибающего момента 
представлены в виде безразмерных величин М”. Они рассчитаны делением размерного изгибающего момента 
на максимальное значение для данного сечения (например для сечения х = 0,3Ё максимальное значение 
Мптах = 4,6 : 10° Нм), найденное при варьировании величины 0 во всем рассматриваемом диапазоне. 
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Малыхина О. И. Аналитическое решение приближенных уравнений движения ракеты-носителя под действием порыва ветра 
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Рис. 4. Зависимость значения безразмерного изгибающего момента М” от продолжительности действия порыва ветра 0”, 
выраженной в долях периода Т! для различных сечений РН: зона 1 (х = 0,0 — 0,2Г,) (а); зона 2 (х = 0,2[ — 0,4Г,) (6); зона 3 
(х = 0,4[ — 0,55Г,) (в); зона 4 (х = 0,55Ё — 0,75, (г); зона 5 (х = 0,75Ё — 0,9Г) (9); зона 6 (х = 0,9[ — Г) (е) 


Итак, продолжительность действия порыва на графиках представлена безразмерной величиной 0”, 
полученной путем деления параметра 0 на период Т, первого тона колебаний РН. Каждая линия графика 
соответствует одному сечению РН. Все сечения РН сгруппированы по характеру функции М”(0^) и приведены 
на различных графиках, а длина РН соответствующим образом разделена на зоны (рис. 5). 
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Рис. 5. Зоны сечений РН 


Из рис. 4 видно, что для первой зоны максимум изгибающего момента достигается уже при значениях 
9 не менее 15 % от периода Т. первого тона колебаний. Для второй и третьей зон — несколько больше 
половины периода Т:. Дальнейшее увеличение значений 0 не влияет на величины максимальных значений 
изгибающего момента. Для зон 4—6 максимум значения изгибающего момента оказывается локальным и 
располагается в зоне значений 09, близких к величине половины периода Т, первого тона колебаний РН. 
Полученные результаты полностью соответствуют итогам конечноэлементного моделирования, проведенного 
ранее”. 

Таким образом, для получения максимальных значений изгибающего момента в сечениях 
рассматриваемой РН необходим динамический анализ поведения конструкции с внешним воздействием в виде 
порыва ветра, продолжительность которого определяется параметром 0, близким по значению к половине 
периода первого тона колебаний РН. 

Обсуждение и заключения. С помощью метода наложений получены аналитические решения, 
описывающие движение системы с одной степенью свободы, на которую влияют сила трения и внешняя сила, 
меняющаяся по трапецеидальному закону. Приведена методика применения полученных аналитических 
решений для систем с большим числом степеней свободы. Показано хорошее совпадение двух видов решений: 

—щ аналитических, для ускорений точек РН, найденных по упрощенной модели; 

— численных, полученных по полной конечноэлементной модели РН. 

Показано, что аналитические решения могут применяться для анализа динамических силовых факторов 
с целью выбора продолжительности порыва ветра, при воздействии которого достигаются максимальные 
нагрузки в сечениях корпуса РН. Аналогично можно анализировать перегрузки, которые достигаются в 
сечениях РН (например, в точках установки систем измерения). 

Кроме того, на основании предложенной методики можно построить полный цикл предварительного 
расчета нагрузок в случае, когда возможно аналитическое представление внешней динамической нагрузки. 
Анализ нагружения на основе аналитических решений весьма экономичен с точки зрения расчетного времени и 
может быть прекрасной альтернативой конечноэлементному моделированию на этапе проектных расчетов, 
когда исследуется большое число вариантов комбинаций внешних нагрузок и конфигураций разрабатываемой 
конструкции. Конечноэлементный анализ подробной модели в этом случае можно использовать как 
уточняющий итоговый расчет. 
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